Rotational intermittency and turbulence induced lift 
experienced by large particles in a turbulent flow. 
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The motion of a large, neutrally buoyant, particle, freely advected by a turbulent flow is 
determined experimentally. We demonstrate that both the translational and angular accelerations 
exhibit very wide probability distributions, a manifestation of intermittency. The orientation of the 
angular velocity with respect to the trajectory, as well as the translational acceleration conditioned 
on the spinning velocity provide evidence of a lift force acting on the particle. 
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The description of a solid object freely advected in a 
fluid requires, in addition to its translational degrees of 
freedom characterizing its position, 3 rotational degrees 
of freedom, specifying its orientation with respect to a 
reference frame. The evolution of its position and of its 
orientation depends, according to Newton's laws, on the 
forces and torques acting on the particle at each instant, 
which result from the interaction between the object and 
the turbulent flow. Their determination raises challeng- 
ing issues. The problem has been solved to a large ex- 
tent for spherical particles of size D much smaller than 
the smallest length scale of the flow, the Kolmogorov 
scale 77 [1, 2]. Because of the small size of the particle, the 
flow around it is locally laminar - see Fig. 1 - therefore 
the equation governing the particle velocity, i;, can be de- 
termined by solving the fluid (Stokes) equations once the 
fluid velocity, n, is known. In the simplest case, the par- 
ticles are subject to the Stokes drag and the added mass 
term [ ] , thus the velocity v can be determined by solving 
a simple differential equation. For D ^ 0, the velocity of 
a neutrally buoyant particle reduces to the fluid velocity, 
n, so the particle behaves as a fluid tracer. This prop- 
erty is crucial for several experimental techniques [ ]. 
Interestingly, the translational and rotational degrees of 
freedom completely decouple in the limit I) ^ [ , ]. 

As suggested by Fig. 1, the case of particles of ra- 
dius larger than r] is conceptually much more difficult. 
Experiments have shown that upon increasing the ratio 
D/t] from 1 to 40, the variance of the particle acceler- 
ation {i.e. of the forces) decreases as {D/r])~'^^^ [ , ]. 
The fluctuations of force are still non-Gaussian up to 
D/f] < 40 ]. A full derivation of the equation of mo- 
tion of a large particle is still not available. The possible 
existence of a lift force when the particle Reynolds num- 
ber Rcp = UrmsD/v is large, coupling rotational and 
translational dynamics, has been the subject of several 
studies [10, 11]. It can be viewed as a generalization of 
the Magnus force, derived in an inviscid, laminar flow: 
^lift = C'liftu;^ X Vrei, where u/ is the rotation of the 
particle in the flow. A lift has been measured in lab- 
oratory experiments, when the flow is steady and lami- 




FIG. 1: Particles in a turbulent flow. Whereas the flow 
around the small particle is smooth, it exhibits significant 
spatial variations around the large particle. 



nar [12, 13]. The flow conditions in these experiments are 
clearly very different from the case of a particle in tur- 
bulence, schematically represented in Fig. 1. The very 
existence of any lift force in these conditions is thus not 
obvious. 

Here, we measure simultaneously the translational 
and rotational motion of a neutrally buoyant spherical 
particle, whose diameter is a fraction of the integral 
size of the turbulence, and report the evidence of a 
lift force. We flrst briefly describe our 6-dimensional 
tracking technique [14-16], then we present our results 
concerning the intermittency of the translational and 
rotational velocity and acceleration. Last, we discuss 
how the translational and rotational degrees of freedom 
of the particle couple, resulting in a lift force. 

The flow and tracking technique are described in de- 
tails in [il]. a von Karman (VK) flow, as widely used 
for turbulence studies [6, 17, 18], is generated in the gap 
between two counter-rotating impellers, see Fig. 2(a). In 
order to be able to perform direct optical measurements, 
the container is built with flat Plexiglas side walls, so 
that the cross section of the vessel has a square shape. 

In the results reported here, the driving disks are ro- 
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FIG. 2: (a) The flow domain in between the impeller has characteristic lengths H — 2R — 20 cm; the disks have radius 
= 10 cm, fitted with straight blades 1 cm in height; (b) It is illuminated by high power LEDs and sequences of 8 bit gray 
scale images are recorded using high-speed cameras; (c) an adequate texture painted on the particle allows the tracking of 
its orientation; (d) Example of particle tracks and orientations (the green and blue arrows mark North-South and East- West 
directions). 



tated at 3 Hz, which corresponds to a power input in 
the system of the order of £ 1.7 W/kg. The integral 
time and length scales are Lint = 3 cm, Tl = 0.3 s, so 
that the dissipative space and time scales are respectively 
77 ~ 30 /im and ~ 1 ms. The flow Reynolds number 
based on the Taylor micro-scale is Rx ~ 300. 

A PolyAmide sphere with density pp = 1.14 g.cm~^ 
and diameter D = 18 mm is made neutrally buoyant 
by adjusting the density of the fluid by addition of 
glycerol to water (the final density mismatch is less 
than Ap / p = 10~^). The spheres are homogeneous, 
i.e. the particle center of mass coincides with its 
geometrical center. Their size is comparable to the 
integral size of the flow, corresponding to D/r] ~ 600 
and D/I/int ^ 0.6. Its motion is recorded at a frame-rate 
equal to 600 Hz, using 2 high-speed video cameras 
(Phantom V12, Vision Research Inc.) positioned at 90 
degrees. The measurement volume lies within 75% of 
the radial and axial distances on either side of the flow 
center, i.e. extends to regions where anisotropy and 
inhomogeneity are known to play a role. The position 
of the sphere is determined using a position tracking 
algorithms, as in [ ]. The time-resolved determination 
of the rotational degrees of freedom [14, 16] is carried 
out by painting a pattern that enables a determination 
of the particle's orientation from the camera images, 
cf. Fig. 2(c). The orientation tracking algorithm then: 
(i) compares the sphere's picture with synthetic images 
and identifies a set of possible orientations; (ii) from 
the set of possible candidates at successive instants, 
a Flow algorithm identifies a likely time series; (iii) a 
post-treatment adjusts remaining ambiguities, using 
the 2 independent views from the 2 cameras. After 
processing, all trajectories with a duration longer than 



0.25 Tl are analyzed - corresponding to 3434 trajecto- 
ries, with duration ranging between 0.25 and 3 T/,. A 
few particular examples are shown in Fig. 2(c), showing 
both the position, color coded for the amplitude of the 
acceleration, together with two unit vectors fixed in the 
reference frame of the particle. 

We find that the behavior of the translational degrees 
of freedom is qualitatively very similar to the results ob- 
tained for much smaller particles [5, 7-9] - the size of the 
particles studied so far did not exceed D < 40/^, whereas 
ours has D ^ 600r]. The particle velocity has a quasi- 
Gaussian distribution; each component has zero mean 
and an rms intensity of fluctuations of about 60 cm/s, 
for impellers rotating at 3 Hz. This is of the order of 
30% of the impeller tip speed, and also of the order of 
the rms velocity fluctuations of tracers in the same flow. 
This leads to the consistent notion that the velocities of 
the particle are of the order of the large scale swirls gen- 
erated in the VK setup. The particle Reynolds number 
Rep = u^msD/v has a mean equal to 1200, with values up 
to 3500; one expects that its Reynolds number relative to 
the local fluid motions is of the same order of magnitude, 
too. 

The particle's acceleration, with zero mean, has an 
rms amplitude of fluctuations equal to 8.5 m/s^. Its 
probability distribution function (PDF) has strongly 
non-Gaussian fluctuations as shown in Fig. 3. In terms 
of velocity increments, the acceleration can be viewed 
as a velocity change over the shortest time scale, while 
the velocity increments over long times scales have the 
same distribution as the velocity itself. Our results thus 
demonstrate that the PDFs of velocity time increments 
evolve from a shape with wide tails [ ] for very small 
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FIG. 3: PDFs of velocity and acceleration, for the linear and 
angular motions. The curves have been shifted vertically for 
clarity. 



time lags to a Gaussian shape at long time intervals, 
a phenomenon often called intermittency. This effect 
is observed here for a very large particle, compared to 
the flow Kolmogorov scale, rj. The strong events in the 
acceleration PDF are, however, less intense than in the 
case of smaller particles. We estimate the normalized 
fourth moment, the flatness to be F = 7 ± 1 for our 
large particle {D/Li^t ~ 0.6). Remarkably, this value is 
identical to the one found in [7] , with a comparable ratio 
D/Lint ^ 1, and generally with the very slow decay of 
the flatness with D observed in [9, 20] for particles much 
smaller compared to Lint {D/Lint 1/50). 

We now turn to the angular velocity of the sphere. 
The three components fluctuate around a zero mean 
value (there is no preferred orientation) and their rms 
amplitude is 12 rad/s. This value is of the order of the 
rotation rate of the driving disks, and also corresponds 
to the rotation that would result from imposing a 
velocity difference of the order of Urms across the particle 
diameter D {u^ms/D 30 rad/s). The PDF of angular 
velocity components are shown in Fig. 3. The distri- 
butions are symmetric, and slightly non- Gaussian (we 
estimate a flatness factor F ~ 4). The rms amplitude 
of the angular acceleration is about 700 rad/s^, again of 
the order of (i^rms/^)^- The statistics of angular acceler- 
ation is strongly non-Gaussian (we estimate F = 7 ± 1). 
Hence, the PDFs of the angular velocity increments 
becomes broader when the time-lag r decreases from 
r ^ Tl to r ~ r^: the angular dynamics is intermittent. 

We now address the question of the coupling between 
the dynamics of rotation and translation. Fig. 4(a) 
reveals a strong alignment between the direction of the 
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FIG. 4: (a) alignment of the particle rotation, u?^, with re- 
spect to the moving Frenet coordinate system; (b) normal 
acceleration conditioned to the component of angular veloc- 
ity parallel to the bi-normal Frenet vector. 



rotation vector lv^ and the vectors defining the trajec- 
tory - the usual Frenet coordinate system (T, A/", B), 
denoting respectively the units vectors along the velocity, 
the curvature, and the direction perpendicular to the 
tangent plane of the trajectory. We find that uj^ is 
aligned perpendicular to AT, so both the acceleration and 
LV^ X V lie in the (T, A/')-plane. While the observation 
that u;^ - N ^ is consistent with a Magnus force, 
the fairly sharp distribution of the direction of oj^ on 
the sphere is unexpected. Fig. 4(a) shows a maximum 
of probability at an angle of the order of 35° with 
respect to T in the (T, S)-plane. The amplitude of the 
acceleration {i.e. force) along the vector AT, conditioned 
on the amplitude of oob = - B is shown in Fig. 4(b). 
The averaged normal acceleration, conditioned on ujb 
increases from 6 to 9 m/s^ (half a standard deviation) 
when the particle rotation varies in the range ±12 rad/s 
{i.e. one standard deviation in rotation speed). 

The main results of this work concern the strong in- 
termittency both for the translational and the angular 
accelerations, and the coupling between the rotational 
and translational degrees of freedom. 

The observed intermittency of the translational and 
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rotational accelerations points to very intense fluctua- 
tions of the force and torque applied to the particle. The 
interaction between the sphere and the flow involves a 
pressure and a viscous terms, resulting from the stress 
tensor: r^j = —pSij + pv{diUj + djUi). Measurements 
of the drag force acting on a particle a large Reynolds 
numbers [zi] demonstrate that pressure is the dominant 
effect in the force: at the particle's Reynolds number in 
our experiment, Rep ~ 10^, the friction force contributes 
of the order of ~ 20% of the drag; this fraction diminishes 
at higher Reynolds numbers. Pressure is expected to be 
mostly coherent at the scale of our particle D ^ so the 
pressure effects act collectively on the object, resulting in 
a force of order pD'^u^^^, which yields the correct order of 
magnitude for the rms of the particle acceleration. The 
strongly non-Gaussian PDF of acceleration fluctuations 
are more surprising, since the force acting on the particle 
results from an averaging over a size comparable to the 
correlation length of the flow. 

The effect of friction, although weak compared to pres- 
sure, is essential to understand the torque. In fact, the 
pressure force, perpendicular to the surface of the sphere, 
does not contribute to the torque applied at the cen- 
ter of the particle. Thus, the rms of the angular ac- 
celeration is (cj^)-^/^ ~ (r^)^/^//, where is the vis- 
cous torque and / the moment of inertia of the particle. 
The magnitude of F^ can be inferred from the estimated 
contribution of the friction force to the drag, Fy/F]j ~ 
vdxU /i/rms ^ ^l/^rms^ whcrc IS the vclocity charac- 
teristic of skin friction velocity, deflned hy u^'^ = vdxU. 
Given the weak value of the moment of inertia tensor, 
one obtains (cj^)-^/^ ~ 10(i^*/i^rms)^('^rms/^)^- With 
{u^/urmsY ^ range of Reynolds numbers of 

the particle motions, this leads to (cj^) ~ ('^?ms/^^)^- 
Which properties of the turbulent flow control the rate of 
rotation of the particle also remains to be elucidated. In 
this respect, the naive vision of many small eddies, com- 
pared to the size of our big particle, implicit in Fig. 1, 
is likely to be an over simpliflcat ion. Small eddies acting 
on the particle in a spatially incoherent manner would 
result in a signiflcantly reduced torque acting on the par- 
ticle. This suggests a much more coherent flow pattern, 
at least at our particle Reynolds numbers, in fact con- 
sistent with the recent numerical results of [22]. We also 
note that while the estimates above provide a qualitative 
explanation of the observed velocities and accelerations 
rms values, understanding the complete statistics of their 
fluctuations remains a challenge. 

Our evidence of the lift effect, i.e. of a strong coupling 
between rotation and translation, rests on conditioning 
the acceleration (force) on the angular velocity of the 
particle, u/' . While this is a very sensible choice, the 
actual torque acting on the particle depends on the local 
structure of the flow, which remains to be determined, 
for instance by local Particle Image (or Tracking) 
Velocimetry around the moving sphere [ ], or numer- 



ically [20, 22, 24]. The sharply peaked distribution of 
u/' in the Frenet basis, see Fig. 4(a), is a very surprising 
result of this work. In this respect, determining experi- 
mentally the variation of the quantities studied here as 
a function of particle size and Reynolds number should 
provide important clues on the interaction between 
turbulent flow and the particle. 
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